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Semiclassical spectra weighted with products of diag- 
onal matrix elements of operators A a , i.e., g aa '(E) ~ 

{n\A a \n){n\A a i\n) /(E — E n ) are obtained by harmonic 
inversion of a cross-correlation signal constructed of classical 
periodic orbits. The method provides highly resolved semi- 
classical spectra even in situations of nearly degenerate states, 
and opens the way to reducing the required signal lengths to 
shorter than the Heisenberg time. This implies a significant 
reduction of the number of orbits required for periodic orbit 
quantization by harmonic inversion. 

PACS numbers: 03.65.Sq, 05.45.+b 



The semiclassical quantization of periodic orbits is of 
fundamental interest for both regular |jj and chaotic [^|J|] 
dynamical systems. In both cases the semiclassical den- 
sity of states is obtained by summing over the periodic 
orbits of the underlying classical system. However, pe- 
riodic orbit sums of this type usually diverge. Various 
techniques have been developed in recent years to over- 
come this problem. Most of these techniques are es- 
pecially designed for hyperbolic (chaotic) systems 
and cannot be applied to systems with regular or mixed 
regular-chaotic dynamics. Recently a new technique, 
based on harmonic inversion of semiclassical signals , 
was shown to be both very powerful and universal for 
the problem of periodic orbit quantization in that it does 
not depend on special properties of the system such as 
ergodicity or the existence of a symbolic dynamics. The 
method only requires the knowledge of periodic orbits 
and their physical quantities up to a certain maximum 
period, which depends on the average local density of 
states. Unfortunately, this method is not free of the gen- 
eral drawback of most semiclassical approaches, which 
suffer from a rapid (exponential for chaotic systems) pro- 
liferation of periodic orbits with their period, which in 
turn requires an enormous number of orbits to be taken 
into account. 

In this paper we propose to implement an extension 
of the filter- diagonalization method, which is a method 
of solving the harmonic inversion problem [9 10 1, to the 
case of time cross-correlation functions ]S Jll| , |l2| , i.e., 
the one-dimensional time signal C(t) is extended to a 
D x D matrix of cross-correlated time signals C aa >(t), 
with a, a' = l,...,D. This method has recently also 



served as a powerful tool for the semiclassical calculation 
of tunneling splittings [ fL3| . The informational content of 
a D x D time signal is increased roughly by a factor of 
D as compared to a 1 x 1 signal. Here we demonstrate 
that the required amount of periodic orbits can be sig- 
nificantly reduced when a cross-correlated periodic orbit 
sum (signal) is inverted instead of a single periodic orbit 
sum. The power of the method is demonstrated for the 
circle billiard, as an example of a completely integrable 
system. 

Consider a quantum Hamiltonian H whose eigenvalues 
are w n and eigenstates \n). We introduce a set of D 
smooth and linearly independent operators [there is no 
need to choose commuting observables] and define a non- 
trace type cross-correlated function fbl 



Jaa'(w) 



IC 



(1) 



where b an and b a ' n are the diagonal matrix elements of 
operators A a and A a i , respectively, i.e. 

b an = (n\A a \n) . (2) 
Note that Eq. [l] can only be written as a trace formula, 
9aa'{w) = tr ^A a G + (w)A a , I 

with the Green Function G + (w) = (w — H + ie)^ 1 , if 
either A a or A a i commutes with H |14|. The weighted 
density of states is given by 



g aa '(w) = lmg aa >(w) 

ir 



(3) 



Until recently the periodic orbit quantization would in- 
volve the use of periodic orbit expressions for functions 
Qaa' (w) similar to that in Eq. The problem one would 
then encounter when trying to extract the semiclassical 
eigenenergies (poles), is that of the analytic continuation 
of g^ a , (w) to the real axis, where the latter diverges. A 
general procedure proposed in Refs. [QJ|] introduced as- 
pects of signal processing into the solution of this prob- 
lem. It was suggested to perform a spectral analysis (har- 
monic inversion) of a "time" (more precisely, action s) 
signal constructed of the same set of periodic orbits. Here 
we extend the harmonic inversion procedure of Refs. [j7]|^] 
to the analysis of a semiclassical cross-correlation signal 
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C aa ,(s) = ^- [ g aa ,(w)e- lsw dw , (4) 

denned as the Fourier transform of g aa '(w). When ap- 
plied to the quantum expression ([[]), this yields 

C a a'{s) — i ^ ' bgnbg'nG n . (5) 
n 

Let us assume that C aa > (s) is given, and the spectral 
parameters (e.g., w n and d aa / n — b an b a i n ) are to be 
extracted. This can be done by solving the conventional 
harmonic inversion problem , which is formulated as 
a nonlinear fit of the signal C(s) by the sum of sinusoidal 
terms, 

n 

with the set of, in general, complex variational parame- 
ters {w ni d n }. Simple information theoretical considera- 
tions ]I(]|,p| then yield an estimate for the required signal 
length, s max ~ Aitg{w), for poles w n < w which can be 
extracted by this method. When a periodic orbit approx- 
imation of the quantum signal C(s) is used, this estimate 
results sometimes in a very unfavorable scaling because 
of a rapid (exponential for chaotic systems) proliferation 
of periodic orbits with increasing period. 

Consider a generalized harmonic inversion problem, 
which assumes that the whole s-dependent D x D sig- 
nal C aa i (s) is adjusted simultaneously to the form of Eq. 
^, with b an and w n being the variational parameters. 
The advantage of using the cross-correlation approach 
|p|, p]JT^1 is based on the simple argument that the to- 
tal amount of independent information contained in the 
D x D signal is D(D + 1) multiplied by the length of the 
signal, while the total number of unknowns (here b an and 
w n ) is (D + l) times the total number of poles w n . There- 
fore the informational content of the D x D signal per 
unknown parameter is increased (compared to the case 
of Eq. |) by a factor of D. [Of course, this scaling holds 
only approximately and for sufficiently small numbers D 
of operators A a chosen.] 

The calculation of a semiclassical approximation to 
C aa ' (s) is significantly simplified for systems with a scal- 
ing property, i.e. where the shape of periodic orbits does 
not depend on the scaling parameter, w, and the classical 
action scales as S po = ws po . For the identity operator 
A\ = I the element Cff (s) is the Fourier transform of 
Gutzwiller's trace formula H|l f° r chaotic systems, and 
of the Berry- Tabor formula [jjj for regular systems, i.e. 
(see Refs. @§) 

Cll(s)=J2A°Hs-s po ) , (7) 

po 

where s po are the periods of the orbits and A po the am- 
plitudes (recurrence strengths) of the periodic orbit con- 
tributions including phase information. For A± = I and 



an arbitrary smooth operator A a the elements Cj^s) 
are obtained from a semiclassical approximation to the 
generalized trace formula tv{G + A a } |E5|S]- The result 
is that the amplitudes A po in ([?]) have to be multiplied 
by the classical average of the observable A a along the 
periodic orbit, 

a a , po = / A a (q(s),p(s))ds , (8) 

Spo Jo 

with A a (q, p) the Wigner transform of the operator A a . 
The problem of finding a semiclassical approximation to 
Eq. ^ for the general case of two arbitrary smooth op- 
erators A a and A a i is investigated in Ref. ]l4|] , where 
numerical evidence is presented that the amplitudes A po 
in (0) have to be multiplied by the product of the classical 
averages, a a , po a a > , po , of the two corresponding classical 
observables, i.e. 

C aa> ( S ) = X! a °<,P° a a' , P oA po 6 (s - S po ) . (9) 
po 

We here adopt the results of Ref. (li) and use Eq. || as 
the starting point for the following application of har- 
monic inversion of cross-correlation functions. Note that 
all quantities in ([)]) are obtained from the classical peri- 
odic orbits. 

The idea of periodic orbit quantization by harmonic 
inversion |7|J^Jll[] is to fit the semiclassical functions 
C^, (s) given in a finite range < s < s max by the func- 
tional form of the quantum expression (j^) . The frequen- 
cies of the harmonic inversion analysis are then identified 
with the semiclassical eigenvalues w n and the amplitudes 
ban with the semiclassical approximations to the diagonal 
matrix elements (niacin). We will show that for a given 
number of periodic orbits the accuracy of semiclassical 
spectra can be significantly improved with the help of 
the cross-correlation approach, or, alternatively, spectra 
with similar accuracy can be obtained from a periodic 
orbit cross-correlation signal with significantly reduced 
signal length. 

Here we only give a qualitative and brief description 
of the method. The details of the numerical procedure 
of solving the harmonic inversion problem (j^) and the 
generalized harmonic inversion problem (|J) are presented 
in Refs. |] [l^] . The idea is to recast the nonlinear fit 
problem as a linear algebraic problem ||. This is done 
by associating the signal C aa >{s) (to be inverted) with 
a time cross-correlation function between an initial state 
Q a and a final state $ Q ' , 

Caa'(s) = <3v|e- iSj ^$ Q ) , (10) 

where the fictitious quantum dynamical system is de- 
scribed by an effective Hamiltonian H c s. The latter is 
defined implicitly by relating its spectrum to the set of 
unknown spectral parameters w n and b an . Diagonaliza- 
tion of H e s would yield the desired w n and b an - This 
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is done by introducing an appropriate basis set in which 
the matrix elements of H c g are available only in terms 
of the known signals C aa '(s). The Hamiltonian H c r is 
assumed to be complex symmetric even in the case of a 
bound system. This makes the harmonic inversion stable 
with respect to "noise" due to the imperfections of the 
semiclassical approximation. The most efficient numeri- 
cal and practical implementation of the harmonic inver- 
sion method with all relevant formulas can be found in 
Refs. @-|l2 . 

We now demonstrate the method of harmonic inver- 
sion of cross-correlated periodic orbit sums for the ex- 
ample of the circle billiard. This is a regular physical 
system, and the periodic orbits and their physical quan- 
tities can be obtained analytically. We choose this sys- 
tem for the sake of simplicity and it will be evident that 
the procedure works equally well with more complex sys- 
tems where periodic orbits have to be searched numeri- 
cally. The nearest neighbor level statistics of integrable 
systems is a Poisson distribution, with a high probability 
for nearly degenerate states, and we will demonstrate the 
power of our new method by fully resolving those nearly 
degenerate states. The exact quantum mechanical eigen- 
values E — h 2 k 2 /2M of the circle billiard are given as 
zeros of Bessel functions J\ m \(kR) — 0, where m is the 
angular momentum quantum number and R, the radius 
of the circle. In the following we choose R = 1. Semi- 
classical eigenvalues can be obtained by an EBK torus 
quantization resulting in the quantization condition |l6| 



kR\J\ — (m/kR) 2 — \m\ arccos - — = ir { n 

rC_ZX 



with m = 0, ±1,±2, ... being the angular momentum 
quantum number and n = 0, 1, 2, . . . the radial quantum 
number. States with angular momentum quantum num- 
ber m/0 are twofold degenerate. 

For billiard systems the scaling parameter is the ab- 
solute value of the wave vector, w = k = |p|/7i, and 
the action is proportional to the length of the orbit, 
S po = hk£ po . The periodic orbits of the circle billiard 
are those orbits for which the angle between two bounces 
is a rational multiple of 2n, i.e., the periods £ po are ob- 
tained from the condition 



£ po = 2m r sin 7 



(12) 



with 7 = irm^/rrir, — 1,2,... the number of turns of 
the orbit around the origin, and m r = 2™^, 2m<£ + 1, . . . 
the number of reflections at the boundary of the circle. 
Periodic orbits with m r ^ 2m^ can be traversed in two 
directions and thus have multiplicity 2. To construct a 
periodic orbit cross-correlation signal C s ^ a , {I) , as defined 
by Eq. 0, we choose three different operators, A% = I 
the identity, A2 — r the distance from the origin, and 
A3 = (L/k) 2 the square of the scaled angular momentum. 



For these operators the classical weights a aiPO (Eq. ||) are 
obtained as 



a l,po — 1 



~ 1 



0-2.no 



a3, P o = cos 7 



cos 7 
tan 7 



arsinh tan 7 



(13) 



The calculation of the weights _4 po in Eq. || depends on 
whether the classical dynamics is regular or chaotic. For 
the circle billiard with regular dynamics we start from 
the Berry- Tabor formula |jj and obtain 



i/4r> 



,3/2 
t-po 



-iff Mpo + f) 



(14) 



where /^ po = 3m r is the Maslov index. Note that 
the formalism is directly applicable to chaotic systems 
with the amplitudes A po in (|9j) computed according to 
Gutzwiller's trace formula 

Once all the ingredients of Eq. || for the circle billiard 
are available, the 3x3 periodic orbit cross-correlation 
signal C^ a , (£) can easily be constructed and inverted 
by the generalized filter-diagonalization method. Results 
obtained from the periodic orbits with maximum length 
Smax = 100 are presented in Fig. |l|. Fig. 0a is part of 
the density of states, g(k), Figs. |l|b and Qc are the den- 
sity of states weighted with the diagonal matrix elements 
of the operators A — r and A = L 2 , respectively. The 
crosses are the results from the harmonic inversion of 
the periodic orbit cross-correlation signals. For compar- 
ison the squares mark the matrix elements obtained by 
exact quantum calculations at positions fc EBK obtained 
from the EBK quantization condition (p"l|). We do not 
compare with the exact zeros of the Bessel functions be- 
cause Eq. |^ is correct only to first order in fi and thus 
the harmonic inversion of C^ a , (s) cannot provide the ex- 
act quantum mechanical eigenvalues. A discussion of the 
semiclassical accuracy can be found in fl6|| , and higher 
order H corrections to the periodic orbit sum are con- 
sidered in [ri}| . However, the perfect agreement between 
the eigenvalues fc HI obtained by harmonic inversion and 
the EBK eigenvalues k EBK is remarkable, and this is even 
true for nearly degenerate states marked by arrows in Fig. 
|l|a. The eigenvalues of some nearly degenerate states are 
presented in Table |[ It is important to emphasize that 
these states with level splittings of, e.g., Ak — 6 x 10~ 4 
cannot be resolved by the originally proposed method of 
periodic orbit quantization by harmonic inversion 
with a periodic orbit signal length s max = 100. To resolve 
the two levels at k « 11.049 (see Table ||) a signal length 
of at least s max w 500 is required if a single periodic or- 
bit function C sc (s) is used instead of a cross-correlation 
function. The method presented in this paper can there- 
fore be used to significantly reduce the required signal 
length and thus the required number of periodic orbits 
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FIG. 1. Density of states weighted with the diagonal ma- 
trix elements of the operators (a) A = I, (b) A — r, (c) 
A = L 2 for the circle billiard with radius R = 1 as functions of 
the wave number k = |p[/7i (in atomic units). Crosses: EBK 
eigenvalues and quantum matrix elements. Squares: Eigen- 
values and matrix elements obtained by harmonic inversion of 
cross-correlated periodic orbit sums. Three nearly degenerate 
states are marked by arrows. 



TABLE I. Nearly degenerate semiclassical states of the 
circle billiard. fc EBK : Results from EBK-quantization. 
fc HI : Eigenvalues obtained by harmonic inversion of 
cross-correlated periodic orbit sums. States are labeled by 
the radial and angular momentum quantum numbers (n,m). 



for periodic orbit quantization by harmonic inversion. As 
such the part of the spectrum shown in Fig. [I] can even 
be resolved, despite the splittings of the nearly degen- 
erate states marked by the arrows, from a short cross- 
correlation signal with s max = 30, which is about the 
Heisenberg period Sh = 2n g(k), i.e. half of the signal 
length required for the harmonic inversion of a 1 x 1 sig- 
nal H . The reduction of the signal length is especially 
important if the periodic orbit parameters are not given 
analytically, as in our example of the circle billiard, but 
must be obtained from a numerical periodic orbit search. 
Note also that the density of periodic orbits increases 
with the period and in chaotic systems the proliferation 
of periodic orbits is exponential. How small can s max get 
as one uses more and more operators in the method? It 
might be that half of the Heisenberg period is a funda- 
mental barrier for bound systems with chaotic dynamics 
in analogy to the Riemann-Siegel formula [|| while for 
regular systems an even further reduction of the signal 
length should in principle be possible. However, further 
investigations are necessary to clarify this point. 

In conclusion, we have introduced a method of periodic 
orbit quantization and calculation of diagonal matrix ele- 
ments based on the construction of a cross-correlated pe- 
riodic orbit sum followed by its harmonic inversion. The 
method is not restricted to bound regular systems but is 
universal and can be applied to open and chaotic systems 
as well. It opens the way to reducing the required signal 
lengths to shorter than the Heisenberg time and there- 
fore significantly reduces the number of orbits required 
for periodic orbit quantization by harmonic inversion. 
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